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1 Introduction 



During many years neutrino physics was a very important branch of elementary particle 
physics. In the last few years the interest to neutrinos particularly increased. This is 
connected first of all with the success of the Super-Kamiokande experiment in which very 
convincing evidence in favour of oscillations of atmospheric neutrinos were obtained. 

It is plausible that tiny neutrino masses and neutrino mixing are connected with the 
new large scale in physics. This scale determines the smallness of neutrino masses with 
respect to the masses of charged leptons and quarks. In such a scenario neutrinos with 
definite masses are truly neutral Majorana particles (quarks and leptons have charges 
and are Dirac particles) It is evident, however, that many new experiments are necessary 
to reveal the real origin of neutrino masses and mixing. 

Experimental neutrino physics is a very difficult and exciting field of research. Now 
it is a time when many new ideas and methods are being proposed. In CERN and other 
laboratories projects of new neutrino experiments are developing. Possibilities of new 
neutrino facility, neutrino factory, are investigated in different laboratories. Thus it is a 
very appropriate time to discuss neutrino physics at the CERN-JINR school. 

In these lectures I will consider different possibilities of neutrino mixing. Then, I will 
discuss in some details neutrino oscillations in vacuum and in matter. In the last part of 
the lectures I will consider the present experimental situation. 

I tried to give in these lectures some important results and details of derivation of 
some results. I hope that lectures will be useful for those who want to study physics 
of massive neutrinos. More results and details can be found in the books Jl|- || and 
reviews ||- [0]- 



Most references to original papers can be found in |17 



2 Neutrino mixing 

According to the Standard Model of electroweak interaction the Lagrangian of the inter- 
action of neutrinos with other particles is given by the Charged Current (CC) and the 
Neutral Current (NC) Lagrangians: 



C? c = -^=£ G W« + h.c, (2.1) 

^ C = -^^J a NC ^- (2.2) 
2 cos Uw 

Here g is the electroweak interaction constant, 6w is the weak (Weinberg) angle and W a 
and Z a are the fields of the and Z° vector bosons. If neutrino masses are equal to 

zero in this case CC and NC interactions conserve electron L e , muon L M and tauon L T 
lepton numbers 

L e = const, = const, L T = const (2.3) 
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Table |2J]. Lepton numbers of neutrinos and charged leptons. 
Lepton numbers of all other particles are equal to zero. 
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The values of the lepton numbers of charged leptons, neutrinos and other particles are 
given in the Table [2.1| . 

According to the neutrino mixing hypothesis masses of neutrinos are different from 
zero and neutrino mass term does not conserve lepton numbers. For the fields of uil that 
enter into CC and NC Lagrangians Q2.1|) and (|2.2j ) we have, in this case, 



VlL = Yl Uli UiL ( 2 - 4 ) 

i 

where is the field of neutrino with mass m 8 and U is the unitary mixing matrix. 



The relation (2.4) leads to violation of lepton numbers due to small neutrino mass 
differences and neutrino mixing. To reveal such effects special experiments (neutrino 
oscillation experiments, neutrinoless double /3-decay experiments and others) are neces- 
sary. We will discuss such experiments later. Now we will consider different possibilities 
of neutrino mixing. 

Let us notice first of all that the relation |2.4j is similar to the analogous relation in 
the quark case. The standard CC current of quarks have the form 

f a c = 2(ull a d' L + c L lJ L + t L lab' L ) (2.5) 

Here 

q=d,s,b q=d,s,b q=d,s,b 

where V is Cabibbo-Kobayashi-Maskawa quark mixing matrix. There can be, however, 
a fundamental difference between mixing of quarks and neutrino mixing. Quarks are 
charged four-component Dirac particles: quarks and antiquarks have different charges. 
For neutrinos with definite masses there are two possibilities: 

1. In case the total lepton number L = L e + L M + L T is conserved, neutrino with definite 
masses are four-component Dirac particles (neutrinos and antineutrinos differ by 
the sign of L) ; 

2. If there are no conserved lepton numbers, neutrinos with definite masses z/j are two- 
component Majorana particles (there are no quantum numbers in this case that can 
allow to distinguish neutrino from antineutrino). 
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The nature of neutrino masses and the character of neutrino mixing is determined by 
the neutrino mass term. 



2.1 Dirac Neutrinos 

If the neutrino mass term is generated by the same standard Higgs mechanism, that is 
responsible for the mass generation of quarks and charged leptons, then for the neutrino 
mass term we have 

£ D = -^^M> IL + h.c. (2.7) 

i,v 

where M D is the complex 3x3 matrix and Vm is the right-handed singlet. In the case 
of mass term ( |2.7[ ) the total Lagrangian is invariant under global gauge invariance 

viL-+e ia is lL , v lR ->e ia is lRl I — > e ia l , (2.8) 

where a is a constant that does not depend on the flavor index I. The invariance under the 
transformation (|2.8|) means that the total lepton number L = L e + L M + L T is conserved 

L = const (2-9) 

Now let us diagonalize the mass term ( |2.7|) . The complex matrix M D can be diagonalized 
by biunitary transformation 

M D = VmU ] , (2.10) 
where V*V = 1, U^U = 1 and = mi5ik, > 0. 



With the help of (|2.10| ), from ( |2.7| ) for the neutrino mass term we obtain the standard 
expression 
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CP = - Y VVR m i( U ^)n u il + h.c. } TUiUiUi (2.11) 



Here 
and 



VIA i=l 



n = v iL + u iR (i = 1,2,3) 



v iR = Ei(V^) a v lR 

For the neutrino mixing we have 



J2u h v lL (2.12) 



Processes in which the total lepton number is conserved, like fi — > e+7 and others, are, 
in principle, allowed in the case of mixing of Dirac massive neutrinos. It can be shown, 
however, that the probabilities of such processes are much smaller than the experimental 
upper bounds. 
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Neutrinoless double /3-decay, 



(A, Z) -> (A,Z + 2) + e~ + e~ , 

due to the conservation of the total lepton number is forbidden in the case of Dirac 
massive neutrinos. 

2.2 Majorana neutrinos 

Neutrino mass terms that are generated in the framework of the models beyond the 
Standard Model, like the Grand Unified SO (10) Model, do not conserve lepton numbers 
L e , L M and L T . Let us build the most general neutrino mass term that does not conserve 
Lei £md L r . 

Neutrino mass term is a linear combination of the products of left-handed and right- 
handed components of neutrino fields. Notice that {vl) C — C{Vl) t is the right-handed 
component and {vr) c = C{Vr) t is the left-handed component. Here C is the charge 
conjugation matrix, that satisfies the relations OyJC -1 = — j a , C T = —C, C^C =10 

The most general Lorentz-invariant neutrino mass term in which flavor neutrino fields 
Vil and right-handed singlet fields vm enter has the following form 



C D - M = --K) c Mfi L + h.c. (2.13) 



Here 




n L = r,j \c with u' L = I Vp,h I and v' R = v^r , (2.14) 




M is complex 6x6 matrix. Taking into account that {vl) c — ~ l/ J i C 1 we have 

C D - M = X -nlC- l Mn L + h.c. . (2.15) 

From this expression it is obvious that there is no global gauge invariance in the case of 
the mass term ( 2.13j ), i.e. that the mass term ( [2.13 ) does not conserve lepton numbers. 



The matrix M is symmetric. In fact, taking into account the commutation properties 
of fermion fields we have 

n T L C- l Mn L = -n T L (C T Y l M T n L = n T L C~ x M T n L . (2.16) 



2 



In fact, L and R components satisfy the relations 

1 + 75 „ 1-75 n 
-2-^ = -a-"* = 

From the first of these relations we have Vl(1 — 7s)/2 = 0. Further, from this last relation we obtain 
[(1 — 75)/2] T F^ = 0. Multiplying this relation by the matrix C from the left and taking into account that 
C^JC" 1 = 7 5 we have [(1 — j^)/2](^l) c — 0. Thus, (vl) C is right-handed component. Analogously we 
can show that (vr) C is left-handed component. 
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From this relation it follows that 

M T = M 

The symmetric 6x6 matrix can be presented in the form 



M 



M L (M D f 
M D Mr 



(2.17) 



where M L = M%, M R = Af£ and M D are 3 x 3 matrices. With the help of (plT]) for the 
mass term (|2.15|) we have 

C D - M = Cf + C D + C^. (2.18) 
Here C D is the Dirac mass term, that we have considered before, and the new terms 



=-^(^a) c M^L + h.c. 



i'i 



(2.19) 



(2.20) 



i', i 



which do not conserve lepton numbers are called left-handed and right-handed Majorana 
mass terms, respectively. The mass term ( 2.13|) is called Dirac-Majorana mass term. 
A symmetrical matrix can be diagonalized with the help of unitary transformation 

M = {U^) T mU ] . 

Here U is unitary matrix and = rrii5ik, m>i > 0. Using the relation (11) we can write 
the mass term (|2.15| ) in the standard form 



1 



1 



£ D -m = — (U*n L ) c mU + n L + h.c. = — vmv 



i 6 



(2.21) 



t=i 



where 



v = U + n L + (U + n L ) 



c 



( v x \ 
\ j 



(2.22) 



Thus the fields (i=l,2...6) are the fields of neutrinos with mass rrii. From ( 2.22 ) it 
follows that the fields Vi satisfy the Majorana condition 



V? = Vi 



(2.23) 



Let us obtain now the relation that connects the left-handed flavor fields un with the 
massive fields v^. From ( ^.22| ) for the left-handed components we have 



n L = Uu L . 



(2.24) 
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From this relation for the flavor field v\l it follows 



^2u H u iL (/ = e,/i,r) 



(2.25) 



i=i 



Thus, in the case of Dirac-Majorana mass term, the flavor fields are linear combinations 
of left-handed components of six massive Majorana fields. From ( |2.25| ) it follows that 
the fields vf R are orthogonal linear combinations of the same massive Majorana fields 



(2.26) 



In the case of Majorana field particles and antiparticles, quanta of the field, are identical. 
In fact, for fermion fields u(x) we have in general case 



u(x) 



1 



1 



(2tt) 3 / 2 J2f 



\p)u r (p)e- ipx + 4(p)C {u r (p)) J 



d?p 



(2.27) 



where c r ip){d r ip)) is the operator of absorption of particle (creation of antiparticle) with 
momentum p and helicity r. If the field u(x) satisfies the Majorana condition ( [2.23j ), then 
we have 

Crip) = d r {p) (2.28) 

Let us stress that it is natural that the neutrinos with definite masses in the case of Dirac- 
Majorana mass term are Majorana neutrinos: in fact there are no conserved quantum 
numbers that could allow us to distinguish particles and antiparticles. 



2.3 The simplest case of one generation (Majorana neutrinos) 

It is instructive to consider in detail the Dirac-Majorana mass term in the simplest case 
of one generation. We have 

^d-m _ _ _ mL (p L y UL _ mV) V R v h - - m R V R (z/r) c + h.c. 



- -in L ) c Mn L + h.c. 



where 



n L = 



M = 



m L m D 
m D m R 



(2.29) 



(2.30) 



Let us assume that the parameters mL,m R and mn are real (the case of CP invari- 
ance). In order to diagonalize the mass term ( 2.29Q let us write the matrix M in the 
form 

M = - TrM + M , (2.31) 
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where Tr M = mi + mo and 



M 



1 



V 



wid -{m R -m L ) ^ 

For the symmetrical real matrix we have 



(2.32) 



Here 



O 



M = OmO T . 



cos d sin $ 
- sin i? cos i? 



is an orthogonal matrix, and m ik = mfiik, where 



m.1,2 = =F- V (m R - m L ) 2 + 4m 2 D 



(2.33) 
(2.34) 

(2.35) 



are eigenvalues of the matrix M . 

From ( p.33| ), (|2.34j ) and (|2.35| ) for the parameters cost? and smfi we easily find the 
following expressions 



cos 2i? 



m R - m L 



a/ ( m R — wil) 2 + 4 m 
For the matrix M from ( |2.33|) and (|2.35|) we have 

M = Om'0 T 

where 



tan 2$ 



2m 



(m fl - m L ) 



(2.36) 



m' 



1.2 -= 2 ( m/? + T V ^ ~ mL ) 2 + 4m D • 



The eigenvalues can be positive or negative. Let us write 

m • = TTliTji , 



(2.37) 



(2.38) 



where = |mj| and 77^ is the sign of the i-eigenvalue. With the help of ( [2.33|) and Q2.38 ) 
we have 

M = (U ] ) T mU ] 

Here 

where y/rj takes the values 1 and i. 

Now using the general formulas ( 2.21|) and ( p . 2 2|) for the mass term we have 



C 



D-M 



2 ^ 

i=l,2 



(2.39) 



Here v% = vf is the field of the Majorana particles with mass mj. The fields vl and (vr) c 
are connected with massive fields by the relation 



vl \ tt Vil 



U( ^ , (2.40) 



v R ) J \ v 2L 

where U = 0(^/rj)* is a 2 x 2 mixing matrix. 
Let us consider now three special cases. 

1. No mixing 

Assume nip = 0. In this case 8 = 0, mi = tul, 1TL2 = rriR and rj = 1 (assuming that 
m L and tur are positive). From fl2.40| ) we have 

v L = v 1L {v R ) c = v 2L . (2.41) 

Thus, if itld = there is no mixing. For the Majorana fields V\ and v 2 we have 

V! = u L + {v L f (2.42) 

v 2 = v R + {v r ) c . (2.43) 

2. Maximal mixing 

Assume uir = rriL, ran ^ 0. From Eq. ( |2.36|) , (|2.37|) and ( p.40| ) we have 

7T 

6* = -, mi,2 = m L ^m D . (2.44) 

(assuming \mn\ < m L ) and 

1 1 . . r 1 1 . . 

Thus if the diagonal elements of the mass matrix M are equal, then we have maximal 
mixing. 

3. See-saw mechanism of neutrino mass generation 
Assume m L = and 

m D < m R (2.46) 
From (|2.35| ) and (|2.37|) we have in this case 

^ Tfi 

nil — — — j m 2 — ttir, ~ (r]i = —1,772 = 1) • (2.47) 

niR m R 

Neglecting terms linear in m^jm^ <C 1, from ( |2.40| ) we have 

v L ^ ~iv 1L , (v R ) c ^ v 2L . (2.48) 
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For the Majorana fields we have 

vx ~ iu L - i(is L ) C \ V2 = vr + (v R ) C ■ (2.49) 

Thus if the condition Q2.46Q is satisfied, in the spectrum of masses of Majorana particles 
there are one light particle with the mass ni\ « ran and one heavy particle with the 
mass mi » nip. The condition m^ = means that the lepton number is violated 
only by the right-handed term —\mRU^{v^) c that is characterized by the large mass 
mR. It is natural to assume that the parameter m^ which characterizes the Dirac term 
—ttidPrPl is of the order of lepton or quark masses. The mass of the light Majorana 
neutrino mi will be in this case much smaller than the mass of lepton or quark. This 
is famous see-saw mechanism. This mechanism connects the smallness of the neutrino 
masses with respect to the masses of other fundamental fermions with violation of the 
lepton numbers at very large scale (usually mp ~ Mqvt — 10 16 GeV. 

In the case of the see-saw for three families in the spectrum of masses of Majorana 
particles there are three light masses m 1 ,m 2 ,m 3 (masses of neutrinos) and three very 
heavy masses Mi, M 2 , M 3 . Masses of neutrinos are connected with the masses of heavy 
Majorana particles by the see-saw relation 

( m}^ 

m,~i-l^< m * (2 = 1,2,3). (2.50) 

where m\ is the mass of lepton or quark in z-family. The see-saw mechanism is a 
plausible explanation of the experimentally observed smallness of neutrino masses. 
Let us stress that if neutrino masses are of the see-saw origin then 

a. neutrinos with definite masses are Majorana particles; 

b. there are three massive neutrinos; 

c. there must be a hierarchy of neutrino masses m,\ <C m 2 <C m%. 



3 Neutrino oscillations 

The most important consequences of the neutrino mixing are so called neutrino oscil- 
lations. Neutrino oscillations were first considered by B. Pontecorvo many years ago in 
1957-58. Only one type of neutrino was known at that time and there was general belief 
that neutrino is a massless two-component particle. B. Pontecorvo draw attention that 
there is no known principle which requires neutrino to be massless (like gauge invari- 
ance for the photon) and that the investigation of neutrino oscillations is a very sensitive 
method to search for effects of small neutrino masses. We will consider here in detail the 
phenomenon of neutrino oscillations. 
Assume that there is neutrino mixing 



} J Uai ViL ■ (3.51) 



where U'U = 1 and z/j is the field of neutrino (Dirac or Majorana) with the mass m^ The 
field v a L in ( |3.51| ) are flavor fields (a = e, //, r) and in general also sterile ones (a = si, ...). 
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Let us assume that neutrino mass differences are small and different neutrino masses 
cannot be resolved in neutrino production and detection processes. 
For the state of neutrino with momentum p we have 



where \ui) is the vector of state of neutrino with momentum p, energy 

Ei = 



2 

mf 



p 2 + m 2 ~ p + (p ^> rrii) . 



(3.52) 



(3.53) 



and (up to the terms mf/p 2 ) helicity is equal to -1. If at the initial time t — the state 
of neutrino is \v a ) at the time t for the neutrino state we have 



(3.54) 



The vector \u a ) is the superposition of the states of all types of neutrino. In fact, from 
Q3.52| ), using unitarity of the mixing matrix, we have 



) = ^2 \ v ^)U a 'i ■ 



(3.55) 



From ( 3.54j) and ( 3.55 ) we have 



where 



A < , Ua (t) = J2 u ^ e ~ lEltu o 



(3.56) 



(3.57) 



is the amplitude of the transition v a — > u a > at the time t. The transition amplitude 
A a '- a {t) has a simple meaning: the term is the amplitude of the transition from the 
state \u a ) to the state |z/,); the term e~ lEit describes the evolution in the state with energy 
E^, the term U a /i is the transition amplitude from the state to the state |^). 

The different gives coherent contribution to the amplitude A u > a;Ua (t) . From ( |3.57D 
it follows that the transitions between different states can take place only if: i) at least 
two neutrino masses are different; ii) the mixing matrix is non-diagonal. In fact, if all 
neutrino masses are equal we have a(t) = e~ tEt UanU^ = e~ tEt 8 a > a . If the mixing 
matrix is diagonal (no mixing), we have A u ' aiUa (t) = e~ lEat 8 a i a . 

Let us numerate neutrino masses in such a way that mi < m<i < ... < m n . For the 
transition probability, from ( |3.57|) , we have the following expression: 



[{e**-** - 1) + 1] U t 
8 aa ' + ^2 Ua>i U ai (e~ tA 



2 -L 

a 2 P — 1 
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where Am-j = mf — m\ and L ~ t is the distance between neutrino source and neutrino 
detector. Thus the neutrino transition probability depends on the ratio ~, the range of 
values of which is determined by the conditions of an experiment. 

It follows from Eq. (|3.58| ) that the transition probability depends in the general case 
on (n — 1) neutrino mass squared differences and parameters that characterize the mixing 
matrix U. The n x n matrix U is characterized by rig = n(n — l)/2 angles. The number 
of phases for Dirac and Majorana cases is different. If neutrino with definite masses Ui 
are Dirac particles the number of phases is equal to = (n — l)(n — 2)/2. If z/j are 

Majorana particles the number of phases is equal to nf 3 = n(n — l)/2. 

Notice that from fl3.5S| ) it follows that transition probability is invariant under the 
transformation 

U ai -> e~ i/3a U ai e ia * (3.59) 

where (3 a and cti are arbitrary real phases. From ( |3.59| ) it follows that the number of 
phases that enter into the transition probability is equal to = {n — l)(n — 2)/2 in 
both Dirac and Majorana cases. We come to the conclusion that additional Majorana 
phases do not enter into the transition probability . Thus, by investigation of neutrino 
oscillations it is impossible to distinguish the case of Dirac neutrinos from the case of 
Majorana neutrinos. 

Let us consider now oscillations of antineutrinos. For the vector of state of antineutrino 
with momentum p from ( |3.51| ) we have 



v, 



u a ) = / J Uai\vi) (Dirac case) (3.60) 

i 

j) = U a i\vi) (Majorana case) (3.61) 



where \Vi) (|^)) is the state of antineutrino (neutrino) with momentum p, energy Ei = 
a/ p 2 + mf ~ p + mf/2p and helicity equal to +1 ( up to mf/p 2 terms). 

In analogy with (|3.57|) for the amplitude of the transition V a —>■ V a i in both Dirac and 
Majorana cases we have 

Av a ,^(t) = j2 u «>i e ~ iEitu ™- ( 3 - 62 ) 



If we compare ( 3.57|) and ( p.62|) we come to the conclusion that 

•Av al -,v a (t) = A^Jt) . (3.63) 

Thus for the transition probabilities we have the following relation 

P(u a -> Va?) = P(z7 Q , -> V a ) . (3.64) 

This relation is the consequence of CPT invariance. If CP invariance in the lepton 
sector takes place then for Dirac neutrinos we have 

U* ai = U ai (3.65) 
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while for Majorana neutrinos, from CP invariance, we have 



U ai Vi = U* ai ; (3.66) 



where rji = ±i is the CP parity of the Majorana neutrino with mass mj. From ( |3.57| ) 
( |3.63| ), fl3.65|) and (|3.66|) it follows that in case of CP invariance we have 



P(u a - v' a ) = P(u a -> V' a ) . (3.67) 

Let us go back to the Eq. ( |3.58| ). It is obvious from ( |3.58| ) that if the conditions of an 
experiment are such that Am^- <C 1 for all i then neutrino oscillations cannot be ob- 
served. To observe neutrino oscillations it is necessary that for at least one neutrino mass 
squared difference the condition Am 2 ^ > 1 is satisfied. We will discuss this condition 
later. 



3.1 Two neutrino oscillations 

Let us consider in details the simplest case of the oscillations between two neutrinos 
v a v a i [a' 7^ a\ a, a' are equal to /jl, e or r, //,...). The index i in Eq. (|3.58|) takes 
values 1 and 2 and for the transition probability we have 

PK - O = \6a> a + U a ,U* a2 (e~ iAm ^ - 1)| 2 (3.68) 



For a' a we have from (p. 68 



P(u a -> u' a ) = P{p a > v a ) = -A a / a (l - cosAm 2 ^) (3.69) 

Here the amplitude of oscillations is equal to 

A a ,., a = 4\U a , 2 \ 2 \U a2 \ 2 (3.70) 

and Am 2 = m 2 — m\. Due to unitarity of the mixing matrix 

\U a2 \ 2 + \U a , 2 \ 2 = I (a'^a) (3.71) 

Let us introduce the mixing angle 9 

\U a2 \ 2 = sin 2 6 \U a , 2 \ 2 = cos 2 9 (3.72) 

Thus the oscillation amplitude A a /. a is equal to 

A a , ;a = sin 2 2# (3.73) 

The survival probabilities P{v a —>■ u a ) and P(zv —>■ u a >) can be obtained from (|3.68| ) or 
from the condition of the conservation of the total probability P(v a — > u a ) + P{v a — ► 
u a i) = 1. We have 

P(u a -> z/ Q ) = P(> a , u a ,) = 1 - i sin 2 20(1 - cos ^^) (3.74) 
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Thus in the case of two neutrinos the transition probabilities are characterized by two 
parameters sin 2 29 and Am 2 . 

Let us notice that in the case of transitions between two neutrinos only moduli of 
the elements of the mixing matrix enter into expressions for the transition probabilities. 
This means that in this case the CP relation ( |3.64j ) is satisfied automatically. Thus, in 
order to observe effects of CP violation in the lepton sector the transitions between three 
neutrinos must take place (this is similar to the quark case: for two families of quarks 
CP is conserved due to unitarity of the mixing matrix). 

We also notice that the expression (|3.69| ) for the transition probability can be written 



in the form 

P(v a -> zv) = i sin 2 29 ( 1 - cos2tt-^- j (3.75) 
2 V A)/ 

where 

L = 4. J- (3.76) 



is the oscillation length. The expression (|3.69|) is written in the units H = c = 1. We can 
write it in the form 



~ sin 2 29 ^1 - cos 2.54Am 2 ^ (3.77) 



P(u a -> u a >) = - 

where Am 2 is neutrino mass squared difference in eV 2 , L is the distance in m (km) and 
E is the neutrino energy in MeV (GeV). For the oscillation length we have 



The Eq. ( 3.69 ) and ( 3.74Q describe periodical transitions (oscillations) between differ- 



ent types of neutrinos due to difference of neutrino masses and to neutrino mixing. The 
transition probability depends periodically on L/E. At the values of L/E at which the 
condition 2.54 Am 2 (L/E) = ix(2n + 1) (n — 0, 1, ...) is satisfied, the transition proba- 
bility is equal to the maximal value sin 2 29. If the condition 2.54 Am 2 (L/E) = 2rm is 
satisfied, the transition probability is equal to zero. 

In order to see neutrino oscillations it is necessary that the parameter Am 2 is large 
enough so that the condition Am 2 (L/E) > 1 is satisfied. This condition allows us to 
estimate the minimal value of the parameter Am 2 that can be revealed in an experi- 
ment on the search for neutrino oscillations. For short and long baseline experiments 
with accelerator (reactor) neutrinos for Am 2 min we have, respectively 10-1 eV 2 , 1CT 2 
- 1CT 3 eV 2 (1CT 1 - 1CT 2 eV 2 , 10~ 2 - 10~ 3 eV 2 ). For atmospheric and solar neutrinos 
for Am 2 min we have 10~ 2 - 10~ 3 eV 2 and 10~ 10 - 1CT 11 eV 2 , respectively. Let us notice 
that in the case of Am 2 (L/E) <C 1, due to averaging over neutrino spectrum and over 
distances between neutrino production and detection points, the term cos Am 2 (L/2p) in 
the transition probability disappears and the averaged transition probabilities are given 
by P(u a — > v^) = \ sin 2 29 and P(v a — > u a ) — 1 — | sin 2 29. 
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3.2 Three neutrino oscillations in the case of neutrino 
mass hierarchy 

The two neutrino transition probabilities (|3.69| ) and ( |3.74j ) are usually used for the anal- 
ysis of experimental data. Let us consider now the case of the transitions between three 
flavor neutrinos. 

General expressions for transition probabilities between three neutrino types are char- 
acterized by 6 parameters and have a rather complicated form. We will consider the case 
of hierarchy of neutrino masses 

rrii WI2 <C WI3 

which corresponds to the oscillations of solar and atmospheric neutrinos (we have in mind 
that Amlx can be relevant for oscillations of solar neutrinos and Am^ can be relevant 
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for oscillations of atmospheric neutrinos; from the analysis of the experimental data it 
follows that Am 2 sol ~ 10~ 5 eV 2 (or 1(T 10 eV 2 ) and Am 2 atm ~ 10~ 3 eV 2 ; see later). We will 
see that transition probabilities have in this case the rather simple two-neutrino form. 

Let us consider neutrino oscillations in experiments for which the largest neutrino 
mass squared difference Am^ is relevant. For such experiments 



Amj 2 - < 1 



2p 



and for the probability of the transition v a 
expression 



(3.79) 

i/ a r, from (|3.58 ) we obtain the following 



P(u n 



6 a ' a + U a , 3 U* 3 le- iAm K- - I 
v a i (a' 7^ a) from (|3.80[) we have 



For the transition probability v 



A 2 L \ 

31 2pJ 



where the amplititude of oscillations is given by 



W a > 3 \ 2 \U, 



q3 I 



(3.80) 



(3.81) 



(3.82) 



Using unitarity of the mixing matrix, for the survival probability we obtain, from 
( EDflp and (|3~82[) , 



P(v n 



where 



-B 



1 — cos Am^ 



2p 



a 3 I 



\u, 



a 3 I 



(3.83) 



(3.84) 



It is natural that Eq. ( ft.81| ) and ( p.82| ) have the same dependence on the parameter 
L/E as the standard two-neutrino formulas and ( |3.74j ): only the largest Am 2 is 



relevant for the oscillations. The oscillation amplitudes A a]a and B a - a depend on the 
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moduli squared of the mixing matrix elements that connect neutrino flavors with the 
heaviest neutrino u 3 . Further, from the unitarity of the mixing matrix it follows that 



\U, 



c-3 I 



\u, 



f t3 1 



+ \u TZ 



1 



(3.85) 



Thus, in three-neutrino case with hierarchy of neutrino masses, the transition probabili- 
ties in experiments for which Am^ is relevant are described by three parameters: Am^ , 
| U e 3 1 2 and | | 2 (remember that in the two neutrino case there are two parameters, Am 2 
and sin 2 26 ). 

Since only moduli of the elements of the mixing matrix enter into transition proba- 
bilities, the relation 

P(u a -> v a ,) = P(V a -> V a >) (3.86) 

holds (as in the two-neutrino case). Thus the violation of the CP-invariance in the lepton 
sector cannot be revealed in the case of three neutrinos with mass hierarchy. Notice that 
the relation 

P{v a — ► u a ) = P(zv -> v a i) , (3.87) 

which takes place in the case of two neutrino oscillations, is not valid in the three-neutrino 
case. 

Let us consider now neutrino oscillations in the case of experiments for which Am^ 
is relevant (Am^ j- > 1). From (|3.57p for the survival probability we obtain in this case 
the following expression 



^|^| 2 e- tAm ^ + |^ 3 | 2 e- lAl 



2 2 =- 



=1,2 



(3i 



Due to averaging over neutrino spectra and source-detector distances, the interference 
term cos Ara^ (L/2p) in Eq. (|3.88|) disappears and for the probability we have 



PilSn 



~A 2 + \u t 



a 3 



(3.89) 



i=l,2 



Further, from the unitarity relation Y^h=i \ U a i\ 2 = 1 we have 



(i-if c ,3r) 2 -2|f/„ 1 | 2 ic/, 



a2 



(3.90) 



i=l,2 



Using ( |3.90| ) we can present the survival probability in the form 



Here 



P(u n 



p(l,2)( 



1 - \U a3 \ 2 ) 2 P^ 2 h 



v a ) + \U, 



a3\ 



1 - - sin 2 26 12 (l - cos 2 Ami — ) 
2 v 21 2p ; 



(3.91) 
(3.92) 



16 



and the angle #i 2 is determined by the relations 



177 I 2 177 I 2 

2 7) u ol • 2 77 u cr2 /n no \ 

cos 0i2 = ^ , TT , 9 , sin 0i2 = ^ , TT , 9 , (3.93) 



Z^i=l,2 \ u ai\ Z-a=l,2 l u <*«l 

The probability P^' 2 )(z/ e — * u e ) has the two-neutrino form and it is characterized by two 
parameters: Am^ and sin 2 2$i 2 . We have derived the expression (|3.92| ) for the case of 



the oscillations in vacuum. Let us notice that similar expression is valid for the case of 
the neutrino transitions in matter. 

The expressions ( |3.81| ), ( J3.83|) and ( |3.92| ) can be used to describe neutrino oscillations 



in atmospheric and long baseline neutrino experiments (LBL) as well as in solar neutrino 
experiments. In the framework of neutrino mass hierarchy, in the probabilities of tran- 
sition of atmospheric (LBL) and solar neutrinos enter different Am 2 (Am 2 x and Am| x , 
respectively) and the only element that connects oscillations of atmospheric (LBL) and 
solar neutrinos is |£/ e 3| 2 . From LBL reactor experiment CHOOZ and Super-Kamiokande 
experiment it follows that this element is small (see later). This means that oscillations of 
atmospheric (LBL) and solar neutrinos are described by different elements of the neutrino 
mixing matrix. 



4 Neutrino in matter 

Up to now we have considered oscillations of neutrinos in vacuum. If there is neutrino 
mixing the effects of the matter can significantly enhance the probability of the transitions 
between different types of neutrinos (MSW effect). We will consider here this effect in 
some details. 

Let consider neutrinos with momentum p. The equation of the motion for a free 
neutrino has the form 

i4f»=Jf„|*M> (4.94) 

Let us develop the state \i/j(t)) over states of neutrinos with definite flavor \u a ) (a = 
e, [A, r). We have 

M*)> =X>aK(t) (4.95) 

a 

where a a (t) is the wave function of neutrino in the flavor representation. From (|4.94j) for 
a a (t) we obtain the equation 

.daJQ = J2(u a \H Q \v a ,)a a ,(t) (4.96) 

a' 

Now we will develop the state \u a ) over the eigenstates \vi) of the free Hamiltonian 



H n : 



FoK) = Ei\vi) , (4.97) 
Ei = Jp 2 + mf ~ p + — L . (4.98) 
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We have: 

i 

If we compare ( |4.99| ) and ( |3.52| ) we find 
Further we have 



2 

{u a \H \u a >) = y](u a \Ui)(Ui\H \Ui)(Ui\u a >) = }UaiTr U ia.> + P S a 



(4.99) 



(4.100) 



(4.101) 



The last term of [4.101J , which is proportional to unit matrix, cannot change the flavor 
state of neutrino. This term can be excluded from the equation of motion by redefining 
the phase of the function a(t). We have: 



.da(t) TT m 2 



Of 



U—rfa(t) 
2p 



(4.102) 



This equation can be easily solved. Let us multiply (|4.102|) by the matrix U* from the 
left. Taking into account unitarity of the mixing matrix we have: 



da' (t) m 2 ... 
i— — = — a(t) 
dt 2p w 



(4.103) 



where a'(t) = Wa(t). The solution of equation ( |4.103| ) has the form 

a '(t) = e-^a'iO) . (4.104) 
For the function a(t) in flavor representation, from (|4 . 1 3|) and ( [4.104! ), we find 

a(t) = Ue'^^aiO) (4.105) 

and for the amplitude of the u a — > u a > transition in vacuum from ( |4.105| ) we obtain the 
expression 

Av ;1 Jt) = J2 U ^ e '^Ki ( 4 -!06) 

i 

which (up to the irrelevant factor e~ vpt ) coincides with ( |3.57| ). 

Let us now introduce the effective Hamiltonian of interaction of flavor neutrino with 
matter. Due to coherent scattering of neutrino in matter, the refraction index of neutrino 
is given by the following classical expression: 



27T 

n{x) = 1 + -=-f(0)p(x) 



(4.107) 
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Here f(0) is the amplitude of elastic neutrino scattering in forward direction, and p(x) is 
the number density of matter (the axis x is the direction of p). The effective interaction 
of neutrinos with matter is determined by the second term of Eq. ( |4. 107 ) : 

2ir 

#/(» = p[n(x) - 1] = — f(0)p(x) (4.108) 

NC scattering of neutrinos on electrons and nucleons (due to the Z-exchange) cannot 
change the flavor state of neutrinos. This is connected with u e — — v T universality of 
NC: the corresponding effective Hamiltonian is proportional to the unit matrix^. 

CC interaction (due to the W-exchange) gives contribution only to the amplitude of 
the elastic v e -e scattering 

v e + e -> v e + e (4.109) 
For the corresponding effective Hamiltonian we have 

Hiix) = -^2v eLl a v eL e la {l - 75 )e + h.c. (4.110) 

The amplitude of process ( |4.109| ) is given by 

i v ^ = ^—G FV (4.111) 
\] 2ix 

and, from (|4.108|) and ( |4.111 ), for the effective Hamiltonian in flavor representation we 
have 

Ei{x) = V2G F p e (x)(3 (4.112) 

where (0) Ve;Ve = 1, while all other elements of the matrix f3 are equal to zero and p e (x) is 
the electron number density at the point x. 

The effective Hamiltonian of the neutrino interaction with matter can be also obtained 
by calculating of the average value of the Hamiltonian ( |4.110| ) in the state which describes 
matter and neutrino with momentum p and negative helicity . Taking into account that 
for non-polarized media 

(mat \e(x)^ a e(x) \ mat) = p e (x)6 a0 , (4.113) 
(mat|e(f)7 a 7 5 e(f)|mat) = 0, (4.114) 

from QOTOD we obtain ((OTJ). 

The evolution equation of neutrino in matter can be written, from ( |4.102|) and ( [012 ), 
in the following form (t = x): 

i da^)_ = (jj^u^ + y/2G FPe {x)p)a(x) (4.115) 
ox 2p 



3 Let us notice that if there are flavour and sterile neutrinos NC interactions with matter must be 
taken into account. 
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Let consider in detail the simplest case of two flavor neutrinos (say, v e and v^). In 
this case we have 

U=( C ° S t Sin tl (4.H6) 
y — sin 17 cos v J 

where 6 is the mixing angle. Further it is convenient to write the Hamiltonian in the 
form 

H = -TrH + H m (4.117) 

where Tr H = ^{m 2 + ?n|) + \p2GFPe- The first term of (|4.117| ), which is proportional 
to the unit matrix, can be omitted. For the Hamiltonian we have then 

l_ ( -Am 2 cos2$ + A(x) Am 2 sin 2tf \ 
{ } Ap V Am 2 sin2tf Am 2 cos2tf - A(x) J \*-*-*-°) 

where Am 2 = m\ — m\ and A(x) = 2y/2Gpp e (x)p. The effect of matter is described by 
the quantity A(x). Notice that this quantity enters only into the diagonal elements of 
the Hamiltonian and has the dimensions of M 2 . 

Let us first consider the case of constant density. In order to solve equation of motion 
we will diagonalize the Hamiltonian. We have: 

R m = l jm E m l jm\ (4.119) 

where E^ 1 is the eigenvalue of the matrix H m and 



It is easy to see that 



E™ 2 = T^fAm 2 cos 26 -A) 2 + (Am 2 sin 26) 2 . (4.121) 

Now, with the help of Eq. ([LTJJ) - Q412TD , for the angle 6 m we have 

Am 2 sin 26 Am 2 cos26»-A 
tan 26 m = - — ; cos 26^ 



Am 2 cos2#-A' y/[ Am 2 cos 26 - A) 2 + (Am 2 sin 2fl) 2 

(4.122) 

The states of flavor neutrinos are given by 

\u e ) = cos6 m \u lm ) +sm6 m \u 2m ); |^) = - sin 6 m \u lm ) + cos 6 m \u 2m ) (4.123) 

where \vi m ) (i — 1,2) are eigenvectors of the Hamiltonian of neutrino in matter and 6 m 
is the mixing angle of neutrino in matter. 
The solution of the evolution equation 

= H m a{x) (4.124) 
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can be now easily found. With the help of Q4.119 ) we have 

.da'(x) 



dx 



E m a'(x) 



where 



a\x) = {U m Ya(x). 



The equation ( |4.125| ) has the following solution: 



a (x) 



e~ iEm ^~^a'(x ) 



(4.125) 
(4.126) 

(4.127) 



where xq is the point where the neutrino was produced. Finally, from fl4.126| ) and ( [4.127] ), 
we have 

a (x) = U m e- iEm{x - X0 \U m ) ] a{x Q ) (4.128) 



The amplitude of the v a —>■ u a > transition in matter turns out to be 

A — \^ Tjrn -iE™{x-x )jj* 

ST-v a i;v a — / J ^a'i u at 



(4.129) 



i=1.2 



and, from Q4.129D and ( |4.120| ), we obtain the following transition probabilities, in full 
analogy with the two-neutrino vacuum case: 



P m (v P 



P m {Ve - Ve) = P™^ 



-sm 2 29 m (l-cosAE m L) 
{l-P m {v e ^v,). 



(4.130) 
(4.131) 



Here AE m = E™ - E™ = ^(Am 2 cos 26 - A) 2 + (Am 2 sin2#) 2 and L = x - x is the 
distance that neutrino passes in matter. 

For the oscillation length of neutrino in matter with constant density we have 



t m 
L 



Air 



P 



^(Am 2 cos 29 - A) 2 + (Am 2 sin 2#) 2 



(4.132) 



The mixing angle and oscillation length in matter can differ significantly from the vacuum 
values. It follows from Q4.122 ) that if the condition^ 



Am 2 cos 29 = A = 2 v / 2G F p f 



P 



(4.133) 



is satisfied, the mixing in matter is maximal (9 m = 7r/4) independently on the value of 
the vacuum mixing angle 9. Notice also that if the condition ( f4.133|) is satisfied, the 
distance between the energy levels of neutrinos in matter is minimal and the oscillation 
length in matter is maximal. We have 

T 

(4.134) 



sin 29 



4 Eq. (4.131) is the condition at which the diagonal elements of the Hamiltonian of neutrino in matter 
vanish. It is evident that in such a case the mixing is maximal. 



21 



where Lq = 4Trp/(Am) is the oscillation length in vacuum. If the distance L in the 
transition probabilities ( |4.131| ) is large (as in the Sun case) the effect of v e — > transitions 
is large even in case of a small vacuum mixing angle 9. The relation ( 4.133 ) is called 
resonance condition. 

The density of electrons in the Sun is not constant. It is maximal in the center of 
the Sun and decreases practically exponentially to its periphery. The consideration of 
the dependence of p e on x allowed to discover possibilities for the large effects of the 
transitions of solar i/ e 's into other states in matter (MSW effect). 

Let us consider the evolution equation when the Hamiltonian depends on the distance 
x that neutrino passes in matter 



,da(x) 
dx 



H m (x)a(x) 



(4.135) 



The Hermitian Hamiltonian H m (x) can be diagonalized by a unitary transformation 

H m (x) = U m (x)E m (x)U m \x) (4.136) 

where U m {x)U m \x) = 1 and E?(x) are eigenvalues of H m (x) . From fl4,135|) and (^4.136|) 
we have 



dt 



where 



a\x) = U m \x)a{x) 
Further, by taking into account that 

U^(x)i^l = + i U-\x) dUm(x] 



dx dx 
we have the following equation for a'(x): 

,da'(x 



dx 



a (X) 



dx 



E m {x) -iU m] (x) 



dU m (x) 
dx 



a (x) 



(4.137) 
(4.138) 

(4.139) 
(4.140) 



In the case p e = const the equation (|4. 140|) coincides with ( [4.125|) . 

Let us now assume that the function p e {x) depends weakly on x and the second term 
in Eq. ( |4.138| ) can be dropped (adiabatic approximation). It is evident that the solution 
of the equation 

■ da i( x ) _ pm,j„/,j (4.141) 



dx 



has the form 



a[(x) 



E?(xM(x) 

c 

E™(x) dx 



a-(zo) (4.142) 

(xq being the initial point). 

It follows from ( |4.141| ) and ( [4.14zi ) that, in the adiabatic approximation, a neutrino 
on the way from the point Xq to the point x remains in the same energy level. From 
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( 4.138|) and ( |4.142| ) we obtain the following solution of the evolution equation in flavor 
representation: 

a(x) = U m (x)e- iJ *o Em{x)dx U m \x )A(X ) . (4.143) 
Moreover the amplitude of v a —>■ v a i transition in adiabatic approximation is given by 



(4.144) 



The latter is similar to the expressions ( |4.106| ) and ( 4.129 ) for the amplitudes of transition 
in vacuum and in matter with p e = const. 
For the case of the two flavor neutrinos 



U m (x) 



costf m (x) sm$ m (x) 
-sm$ m (x) cos$ m (x) 



and tan2# m (x) and cos26 m (x) are given by Eq. ( |4.122| ) in which 

A(x) = 2V2G FPe (x)p 



(4.145) 



(4.146) 



The eigenvalues of the Hamiltonian H m (x) are given by Eq. ( f4.121| ). From (|4.145|) we 
have 

/ d9 m {x) \ 

dU m (x) _ | 

d9 m {x) 



U m Hx) 



dx 



\ 



dx 



dx 




(4.147) 



and the exact equation (^4.140 ) takes the form 



.d_ 

dx \ a'-, 



( 



dx 



J 



(4.148) 



V dx 

The Hamiltonian H m in the right-hand side of this equation can be written in the form 



H m = -(Er + E™) + 



\ 



QQm 

dx 



dx 
hAE m 



(4.149) 



where AE m = E™ — E™. As we stressed several times, the term of the Hamiltonian 
which is proportional to the unit matrix is not important for flavor evolution. 

From Eq. (|4. 149|) it follows that adiabatic approximation is valid if the condition 



d9 r ' 



dx 



< -AE m 



(4.150) 



is satisfied. With the help of ( (4.122j ) it is easy to show that ( f4.150| ) can be written in the 
form 



m 2 sin 29 



dp e 
dx 



< [(Am 2 cos 29 - A) 2 + (Am 2 sin 2^) 2 ] 



21 3/2 



(4.151) 
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If the resonance condition 

Am 2 cos 26 = A(x R ) (4.152) 

is satisfied at the point x = Xr, the condition of validity of the adiabatic approximation 
can be written in the form 



2pcos 29 \-^lnp e (x R )\ 
Am 2 sin 2 26 



<1. (4.153) 



From Eq. Q4.144 ) we obtain the following probability for the v a — > v a > transition in 



the adiabatic approximation: 

P(u a - vj) = \U2i(x)\ 2 \U™(x )\ 2 + (4.154) 

i 

+2i? e ^^(x)^* e - i ^^-^^C/r(xo)^(^o). 

i<k 

For solar neutrinos the second term in the r.h.s. of this expression disappears due to 
averaging over the energy and the region in which neutrinos are produced. Hence for the 
averaged transition probability we have 

P {Va _ Val) = \U^(x)\ 2 \U™(x )\ 2 (4.155) 

i 

Thus, in the adiabatic approximation, the averaged transition probability is determined 
by the elements of the mixing matrix in matter at the initial and final points. For the 
case of two neutrino flavors we have the following simple expression for the v e survival 
probability 

P(u e -> u e ) = cos 2 6 m (x) cos 2 6 m (x ) + sin 2 6 m (x) sin 2 6 m (x ) 

= - (1 + cos2^ m (s) cos2^ m (s )) (4.156) 



From Eq. ( |4. 156| ) it is easy to see that if the neutrino passes the point x = xr where the 



resonance condition is satisfied, a large effect of disappearance of v e will be observed. In 
fact, the condition ( |4.152| ) is fulfilled if cos 29 > (neutrino masses are labelled in such 



a way that Am 2 > 0). At the production point Xq the density is larger than at point 
x R and A(xo) > Am 2 cos29. From ( |4.122| ) it follows than cos2#(x ) < 0. Thus, if the 
resonance condition is fulfilled, we see from Eq. ( |4.156| ) that P{y e — > u e ) < |. If the 
condition 

A{x ) > Am 2 (4.157) 

is satisfied for neutrinos produced in the center of the Sun, then cos26' m (a;o) — —1 and, 
for neutrinos passing through the Sun, the survival probability is equal to: 

P(z/ e ^z/ e )~i(l-cos2#) (4.158) 



24 



It is obvious from this expression that the u e survival probability at small 9 is close to 
zero: all u e 's are transformed into z/^'s. 

Let us consider evolution of neutrino states in such a case. From Eq. ( 4.122 ) it follows 
that, at the production point, 9 m (xo) — tt/2. From (|4.123|) we have then 

We) - \v2m) ; k/i) = ~Wlm) (x = Xq) (4.159) 

Thus at the production point flavor states are states with definite energy. In the adiabatic 
approximation there are no transitions between energy levels. In the final point p e = 
and at small 9 we have 

W2) ^ K), Wi) ^ We) (x = x ) (4.160) 

Thus, all z/ e 's transfer to u^s. The resonance condition ( |4.152|) was written in units 
h = c = 1. We can rewrite it in the following form 

Am 2 cos 29 ~ 0.7 • 10~ 7 EpeV 2 

where p is the density of matter in g- cm -3 and E is the neutrino energy in MeV. In 
the central region of the Sun p ~ 10 2 g • cm -3 and the energy of the solar neutrinos is 
~ lMeV. Thus the resonance condition is satisfied at Am 2 ~ 10~ 5 eV 2 . 

The expression ( |4.155|) gives the averaged survival probability in the adiabatic ap- 
proximation. In the general case we have 

p {Ua ^ Val) = \U^x)\ 2 P lk \U™{x )\ 2 (4.161) 

where is the probabilty of transition from the state with energy E™ to the state with 
energy E™. Let us consider the simplest case of the transition between two types of 
neutrinos. From the conservation of the total probability we have 

Pu = l-P 2 i, P 22 = l-Pi2, ^12 = ^21 (4.162) 

Thus in the case of two neutrinos all transition probabilities are expressed through 
P\2- With the help of ( |4.145| ), ( f4.161[ ) and ( [4.162[ ), for the v e survival probability we have: 

P(u e -f u e ) = I + f l -P 12 J cos2# m (a:)cos2# m (a;o) (4.163) 



2 V 2 

In the literature there exist different approximate expressions for the transition prob- 
ability P 12 . In the Landau-Zenner approximation, based on the assumption that the 
transition occurs mainly in the resonance region, 

P 12 = e -f^ (4.164) 

where 

jAE m _ Am 2 sin 2 29 

lR ~ \d9^/dx\ ~ 2pcos29\f x \np e {x R )\ ' (4 ' 165) 
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In the above equation F = 1 for linear density and F = 1 — tan 2 9 for exponential density. 
The adiabatic approximation is valid if 7^ 3> 1 [see ( |4.150| )1. In this case Pyi — 0. 

This concludes the considerations on the phenomenological theory of neutrino mixing 
and on the theory of neutrino oscillations in vacuum and in matter. We will start now 
the discussion of experimental data. There are three methods to search for the effects of 
neutrino masses and mixing: 

I. The precise measurement of the high energy part of /3-spectrum; 

II. The search for neutrinoless double /3-decay; 

III. The investigation of neutrino oscillations . 

We shall discuss now the results which have been obtained in some of the most recent 
experiments. 

5 Search for effects of neutrino mass 
in experiments on the measurement 
of the /3-spectrum of 3 H 

We will discuss here briefly the results of searching for effects of neutrino masses in 
experiments on the measurement of the high-energy part of the /3-spectrum in the decay 

3 H^ 3 He + e- + z^ (5.166) 



The process Q5.166| ) is a superallowed /3-decay: the nuclear matrix element is constant 



and the /3-spectrum is determined by the phase-space factor and the Coulomb interaction 
of the final e~ and 3 He. For the /^-spectrum we have 

^ = CpE(Q - T)y/{Q-T)*-mlF{E) (5.167) 

Here p is electron momentum, E = m e + T is the total electron energy, Q = m3H — 
m3He — fn e — 18.6 keV is the energy release, C = const and F(E) is the Fermi function, 
which describes the Coulomb interaction of the final particles. In the Eq. ([5.167 ) the 



term (Q — T) is the neutrino energy (the recoil energy of 3 He can be neglected) and the 
neutrino mass enters through the neutrino momentum p u = a/(Q — T) 2 — m 2 . Notice 
that in the derivation of Eq. ( |5.167| ) the simplest assumption was done that v e is the 



particle with mass m v . 

The Kurie function is then determined as follows 



UN 1 



K{T) = ^^FjW) = >fCyl{Q - TW{Q-Tf-rnl (5.168) 

If m u = 0, the Kurie function is the stright line K(T) = \/~C(Q — T) and T max = 0. If 
m v ^ then T max = Q — m u and at small m v the Kurie function deviates from the stright 
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Table 



5.2. 



Neutrino mass from 3 H experiments. 



Experiment 


ml 




m v 




Troitsk 


-1.0 ±3.0 ±2.0 


eV 2 


< 2.5 


eV 


Mainz 


-0.1 ±3.8 ± 1.8 


eV 2 


< 2.8 


eV 



line in the region close to the maximum allowed energy. Thus, if m u ^ in the end point 
part of the spectrum a deficit of observed events must be measured (with respect to the 
number of events expected at m v = 0). 

In experiments on the search for effects of neutrino mass by 3 H-method no positive 
indications in favour of m v ^ were found. In these experiments some anomalies were 
observed. First, practically in all experiments the best-fit values of m 2 are negative. 
This means that instead of a deficit of events, some excess is observed. Second, in the 
Troitsk experiment a peak in electron spectrum is observed at the distance of a few eV 
from the end. The position of the peak is changed periodically with time. There are 
no doubts that new, more precise experiments are necessary. The results of two running 
experiments are presented in Table |5.2| . 



6 Neutrinoless double /3-decay 

The decay 

(A,Z) -> (A,Z + 2) + e~ + e~ (6.169) 

is possible only if the total lepton number L is not conserved, i.e. if neutrinos with definite 
masses are Majorana particles. There are many experiments in which neutrinoless double 
/5-decay ((/3/3) „-decay) of 76 Ge, 136 Xe, 130 Te, 82 Se, 100 Mo and other even-even nuclei is 
searched for. 

Let consider the process ( |6.169| ) in the framework of neutrino mixing. The standard 
CC Hamiltonian of the weak interaction has the form 

Hj = 9z 2 e Ll a u eLja + h.c. (6.170) 
Here j a is the weak hadronic current and 

VeL = UeML (6.171) 

where Vi is the Majorana neutrino field with mass rrii. 

The (PP)ou decay is a process of second order in Gf with an intermediate virtual 
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neutrino. Neutrino masses and mixing enter into the neutrino propagator^ 



(2?r) 4 J - 771 

Taking into account that 

2 (p + to.) g =mi , (6.173) 

we come to the conclusion that the matrix element of (/?/9)oy-decay is proportional tcf] 

<m>=^U 2 ei m i (6.174) 



From ( |6.173|) it is evident that the proportionality of the matrix element of (/3/3)o^-decay 



to < m > is due to the fact that the standard CC interaction is the left-handed one. 
If neutrino masses are equal to zero (/3/3) 0i ,-decay is forbidden (conservation of helicity). 
Notice that, if there is some small admixture of right-handed currents in the interaction 
Hamiltonian, the L — R interference gives a contribution proportional to the j> term in the 
neutrino propagator. Other mechanisms of (/3/3)o-decay are also possible (SUSY with 
violation of R-parity ect.). 

In the experiments on the search for ( / 9 / 9) 0y -decay very strong bounds on the life-time 
of this process were obtained. The results of some of the latest experiments are presented 
in Table |H| From these data upper bounds for | < 777 > | can be obtained. The upper 
bounds depend on the values of the nuclear matrix elements, the calculation of which is 
a complicated problem. From 76 Ge data it follows 

I < m > I < (0.5 - 1) eV (6.175) 

In the future experiments on the search for (/3/3)o„-decay (Heidelberg-Moscow, NEMO, 
CUORE and others) the sensitivity | < m > | < 0.1 eV is planned to be achieved. 



7 Neutrino oscillation experiments 

We will discuss now the existing experimental data on the search for neutrino oscillations. 
There exist at present convincing evidences in favour of neutrino oscillations, which were 
obtained in atmospheric neutrino experiments and first of all in the Super-Kamiokande 
experiment. Strong indications in favour of neutrino masses and mixing were obtained in 

5 We have used the relation vj = —ViC that follows from the Majorana condition — Cvf = z^. 
It is obvious that in the case of Dirac neutrino the propagator is equal to zero. 

6 The term in denominator is small with respect to characteristic p in nuclei (~ 10 MeV) and can 
be neglected. 
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Table Lower bounds of the life-time Ti/ 2 of (/3/3)o^-decay 



Experiment 


Element 


Lower bound of T\m 


Heidelberg- Moscow 


76 Ge 


> 1.6 x 10 25 y 


Caltech-PSI-Neuchatel 


136 Xe 


> 4.4 x 10 23 y 


Milano 


iso Te 


> 7.7 x 10 22 y 



all solar neutrino experiments. Finally, some indications in favour of — > v e transitions 
were obtained in the LSND accelerator experiment. In many reactor and accelerator short 
baseline experiments and in the reactor long baseline experiments CHOOZ no indication 
in favour of neutrino oscillations was found. We will start with the discussion of the 
results of solar neutrino experiments. 

7.1 Solar neutrinos 

The energy of the Sun is generated in the reactions of the thermonuclear pp and CNO 
cycles. The main pp-cycle is illustrated in Fig.[7.1|. The energy of the sun is produced in 
the transition 

4p + 2e" -f 4 He + 2z/ e , (7.176) 

If we assume that solar v e 's do not transfer into other neutrino types [P[y e —* v e ) = 1) 
we can obtain a relation between the luminosity of the Sun, L and the flux of solar 
neutrinos. Let us consider neutrino with energy E. From ( |7. 1 76| ) it follows that 

\{Q~m (7.177) 

is the luminous energy corresponding to the emission of one neutrino. Here 

Q = Am p + 2m e - m iRc ~ 26.7MeV (7.178) 



is the energy release in the transition ( 7.176 ). If we multiply ( [7.1771 ) by the total flux of 



solar z/ e 's from different reactions and integrate over the neutrino energy E we will obtain 
the flux of luminous energy from the Sun 

\ {{Q-2E)Y,h{E)dE = ^j- 2 . (7.179) 

Here L Q ~ 3.86 • 10 33 erg/s is the luminosity of the Sun, R is the Sun-Earth distance and 
If{E) is the flux of neutrinos from the source i (i = pp, ...). Notice that in the derivation 
of the relation (|7.179|) we have assumed that the Sun is in a stationary state. 

The luminosity relation (|7.179 ) is solar model independent constraint on the solar 



neutrino fluxes. The flux I%{E) can be written in the form 

h(E) = Xi(E)$i (7.180) 
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(pp) 



p + p -> 2 H + e+ + v e 



p + e + p^ 2 H + u e 



(pep) 



99.6% 



0.4% 



85% 



3 He + 3 He^ 4 He + 2p 



2 H + p -> 3 He + 7 




lo% 



3 He + 4 He -> 7 Be + 7 



2 x 10" 5 % 



3 He + p -> 4 He + e+ + y e 



(hep) 



99.87% 



0.13% 



( 7 Be) 



7 Be + e~ -> 7 Li + y e 



7 Li + p^2 4 He 



7 Be + p -> 8 B + 7 



8 B -5- 8 Be* + e+ + J/„ 



( 8 B) 



8 Be* -y 2 4 He 



Figure The cycle (the figure is taken from 
ref. 
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Table 



VI- 



Main sources of solar i/ e 's. 



Reaction 


Maximal energy 


Standard Solar Model flux (cm 2 s 1 ) 


pp — > d e + v e 


^ 0.42 MeV 


6.0 x 10 10 


e~ 7 Be — * v e 7 Li 


0.86 MeV 


4.9 x 10 9 


8 B -> 8 Be e+ z/ e 


^ 15 MeV 


5.0 x 10 6 



where $j is the total flux and the function Xi(E) describes the form of the spectrum 
(J Xi(E)dE = 1). The functions Xi(E) are known functions, determined by the weak 
interaction. The luminosity relation (|7.179|) can be written in the form 



i - 2 ir i= ^p (7l8i; 



where E { = J EXi(E)dE is the average energy of neutrinos from the source i. The main 



sources of solar neutrinos are listed in Table 1774 



As it is seen from the Table, the main source of solar neutrinos is the reaction p + 
p — > d + e + + v e . This reaction is the source of low energy neutrinos. The source of 
monochromatic medium energy neutrinos is the process 

e~ + 7 Be -> u e + 7 Li. (7.182) 

The reaction 8 B — > 8 Be + e + + v e is the source of the rare high energy neutrinos. The 
results of solar neutrino experiments are presented in Table [7.5| . 

Homestake, GALLEX and SAGE are radiochemical experiments. In the Kamiokande 
and the Super-Kamiokande experiments recoil electrons (angle and energy) in the elastic 
neutrino-electron scattering are detected. In these experiments the direction of neutrinos 
is determined and it is confirmed that the detected events are from solar neutrinos. 

In the Homestake experiment, because of high threshold [Eth — 0.81 MeV) mainly 8 B 
neutrinos are detected: ~ 77% of events are due to 8 B neutrinos and ~ 15% of events are 
due to 7 Be neutrinos. In GALLEX and SAGE experiments [Eth — 0.23 MeV) neutrinos 
from all reactions are detected: ~ 54% of events are due to pp neutrinos, ~ 27% of 
events are due to 7 Be neutrinos and ~ 10% of events are due to 8 B neutrinos. In the 
Kamiokande and Super-Kamiokande experiments due to the high threshold [Eth — 7 MeV 
for Kamiokande and E t h = 5.5 MeV for the Super-Kamiokande) only high energy 8 B 
neutrinos are detected. 



The results of the solar neutrino experiments are presented in Table |775| As it is seen 
from the Table, the detected event rates in all solar neutrino experiments are significantly 
smaller than the predicted one.[] The most natural explanation of the data of solar 

7 Notice that in the framework of neutrino oscillations the possibility of deficit of solar v^s was 
discussed by B. Pontecorvo in 1968 before the results of the Homestake experiment were obtained. 
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Table |7.5| . Results of solar neutrino experiments. 
[ 1 SNU = 1(T 36 events/ (atoms • sec) ] 



Experiment 


Observed rate 


Expected rate 


Homestake 






v e 37 C1 -> e" 37 Ar 


2.56 ±0.16 ±0.16 SNU 


7.7 ± 1.2 SNU 


Eth = 0.81 MeV 






GALLEX 






z/ e 71 Ga^ e- 71 Ge 


77.5 ± 6.2lt? SNU 


129 ± 8 SNU 


E th = 0.23 MeV 






SAGE 






z/ e 71 Ga^ e" 71 Ge 


66.6 ± If;? + 3 ;^ SNU 


— • — 


E tft = 0.23 MeV 






Kamiokande 






z/e — > z/e 


(2.80 ± 0.19 ± 0.33) 10 6 cm- 2 s- 1 


(5.15^;°°) lOWV 1 


Eth = 7.0 MeV 






Super — Kamiokande 






z/e — > z/e 


(2.44 ± 0.05ig;g?) l^cm-V 1 




E th = 5.5 MeV 







neutrino experiments can be obtained in the framework of neutrino mixing. In fact, 
if neutrinos are massive and mixed, solar z/ e 's on the way to the earth can transfer 
into neutrinos of the other types that are not detected in the radiochemical Homestake, 
GALLEX and SAGE experiments. In Kamiokande and Super-Kamiokande experiments 
all flavor neutrinos u e , and u T are detected. However, the cross section of (u T ) —e 
scattering is about six times smaller than the cross section of v e — e scattering. 

All existing solar neutrino data can be explained if we assume that solar neutrino fluxes 
are given by the Standard Solar Model (SSM) and that there are transitions between two 
neutrino types determined by the two parameters: mass squared difference Am 2 and 
mixing parameter sin 2 26. We will present the results of such analysis of the data later 
on. 

Now we will make some remarks about a model independent analysis of the data. 
First of all from the luminosity relation ( |7.17£ ) for the total flux of solar neutrinos we 
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have the following lower bound 



i 

Furthermore, for the counting rate in the gallium experiments we have 

Qca = [ a(E)J2h(E)dE = J2°i®i = (76±2)SNU (7.184) 

JE th i 

By comparing this lower bound with the results of the GALLEX and SAGE experiments 



(see Table [7.5|) we come to the conclusion that there is no contradiction between experi- 
mental data and luminosity constraint if we assume that there are no transitions of solar 
neutrinos into other states {P{y e — > u e ) = 1). 

It is possible, however, to show in a model independent way that the results of different 
solar neutrino experiments are not compatible if we assume P[y e — > v e ) = 1. In fact, let 
us compare the results of the Homestake and the Super-Kamiokande experiments. We 
will consider the total neutrino fluxes $j as free parameters. From the results of Super- 
Kamiokande experiment we can determine the flux of 8 B neutrinos, $8£ (see Table f7~5"D . If 
we calculate now the contribution of 8 B neutrinos into the counting rate of the Homestake 
experiment we get 

Qc B i = (2.78 ± 0.27) SNU (7.185) 

The difference between measured event rate and Q^f gives the contribution to the Chlo- 
rine event rate of 7 Be and other neutrinos. We have 

Q 7 c f + - = - Q 8 ^ = (-0.22 ± 0.35) SNU (7.186) 

All existing solar models predict much larger contribution of 7 Be neutrinos to the 
Chlorine event rate: 

Qcf(SSM) = (1.15 ± 0.1)SNU (7.187) 

The large suppression of the flux of 7 Be neutrinos (together with the observation of 8 B 
neutrinos) is the problem for any solar model. The 8 B nuclei are produced in the reaction 
p + 7 Be — > 8 B + 7 and in order to observe neutrinos from 8 B decay enough 7 Be nuclei 
must exist in the Sun interior. We can come to the same conclusion about the suppression 
of the flux of 7 Be neutrinos if we compare the results of Gallium and Super-Kamiokande 
experiments. 

All existing solar neutrino data can be described if there are oscillation between two 
neutrino flavors, the neutrino fluxes being given by the SSM values. If we assume that 
the oscillation parameters Am 2 and sin 2 26 are in the region in which matter MSW effect 
can be important, then from the fit of the data two allowed regions of the oscillation 
parameters can be obtained. For the best fit values it was found 



Am 2 = 5 ■ 10~ 6 eV 2 sin 2 26 = 5 ■ 10" 3 (SMA) (7. 
Am 2 = 2 ■ 10" 5 eV 2 sin 2 26 = 0.76 (LMA) (7.189) 
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The data can be also described if we assume that the oscillation parameters are in the 
region in which matter effects can be neglected (the case of vacuum oscillations). For the 
best fit values it was found in this case 

Am 2 = 4.3 • 10~ 10 eV 2 sin 2 2fl = 0.79 (VO) . (7.190) 

In the Super-Kamiokande experiment during 825 days 11240 solar neutrino events 
were observed. Such a large statistics allows the Super-Kamiokande collaboration to 
measure the energy spectrum of the recoil electrons and day/night asymmetry. No sig- 
nificant deviation from the expected spectrum was observed (may be with the exception 
of the high energy part of the spectrum). For the day/night asymmetry the following 
value was obtained 

1 fN-D\ , , 

0.065 ±0.031 ±0.013 (7.191) 



2 \N + D / 

From the analysis of the latest Super-Kamiokande data the following best-fit values of 
the oscillation parameters were found: 

Am 2 = 5 ■ 10" 6 eV 2 sin 2 29 = 5 ■ 10~ 3 (SMA) (7.192) 

Am 2 = 3.2 • 10" 5 eV 2 sin 2 29 = 0.8 (LMA) (7.193) 

Am 2 = 4.3 • 10" 10 eV 2 sin 2 29 = 0.79 (VO) (7.194) 

These values are compatible with the ones in Eq. ( |7.188|) , ( [7.1 89|) and (|7.190|) , which were 
found from the analysis of the event rates measured in all solar neutrino experiments. 

The new solar neutrino experiment SNO started recently in Canada. The target in 
this experiment is heavy water (1 kton of D2O) and Cerenkov light is detected by ~ 10 
photomultipliers. Neutrinos will be detected through the observation of the CC reaction 

v e + d^e~+p + p (7.195) 



as well as of the NC reaction 
and v — e elastic scattering 



u+d^u+n+p (7.196) 
v + e -> v + e (7.197) 



The detection of neutrinos via the CC process ( |7.195| ) will allow to measure the spectrum 
of v e on the Earth. The detection of neutrinos via the NC process ( 7.196| ) (neutrons will 
be detected) will allow to determine the total flux of flavor neutrinos z/ e , u^, v T . From the 
comparison of NC and CC event rates model independent conclusions on the transition 
of solar i/ e 's into other flavor states can be made. 

Next solar neutrino experiment will be BOREXINO. In this experiment 300 tons 
of liquid scintillator of very high purity will be used. Solar neutrinos will be detected 
through the observation of the recoil electrons in the process 

v + e -> v + e. (7.198) 

The energy threshold in the BOREXINO experiment will be very low, about 250 keV. 
That will allow to detect the monoenergetic 7 Be neutrinos. If vacuum oscillations are the 
origin of the solar neutrino problem, a seasonal variation of the 7 Be neutrino signal (due 
to excentricity of the Earth orbit) will be observed. 
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7.2 Atmospheric neutrinos 

Atmospheric neutrinos are produced mainly in the decays of pions and muons 

7i —> \i + 11 — »• e + v e + Uft (7.199) 

pions being produced in the interaction of cosmic rays in the Earth atmosphere. Notice 
that in the existing detectors neutrino and antineutrino events cannot be distinguished. 
At small energies, < 1 GeV, the ratio of fluxes of u^'s and u e 's from the chain ( 7.199]) is 



equal to two. At the higher energies this ratio is larger than two (not all muons decay 
in the atmosphere) but it can be predicted with accuracy better than 5% (the absolute 
fluxes of muon and electron neutrinos are predicted presently with accuracy not better 
than 20 - 25%). This is the reason why the results of the measurements of total fluxes 
of atmospheric neutrinos are presented in the form of a double ratio 

R = wSt (7 - 200) 

where (iV^/iVjdata is the ratio of the total number of observed muon and electron events 
and {N^/ N £ )mc is the ratio predicted from Monte Carlo simulations. 

We will discuss the results of the Super-Kamiokande experiment. A large water 
Cerenkov detector is used in this experiment. The detector consists of two parts: the 
inner one of 50 kton (22.5 kton fiducial volume) is covered with 11146 photomultipliers 
and the outer part, 2.75 m thick, is covered with 1885 photomultipliers. The electrons and 
muons are detected through the observation of the Cerenkov radiation. The efficiency 
of particle identification is larger than 98%. The observed events are divided in fully 
contained events (FC) for which Cerenkov light is deposited in the inner detector and 
partially contained events (PC) in which the muon track deposits part of its Cerenkov 
radiation in the outer detector. FC events are further divided into sub-GeV events (E vis < 
1.33 GeV) and multi-GeV events E v i s > 1.33 GeV). In the Super-Kamiokande experiment 
for sub-GeV events and multi-GeV events (FC and PC) the following values of the double 
ratio R were obtained, respectively (848.3 days): 



R = 0.680^| ± 0.053 



R = 0.6781^ ± 0.080 



(7.201) 



These values are in agreement with the values of R obtained in other water Cerenkov 
experiments (Kamiokande and 1MB) and in the Soudan2 experiment in which the detector 
is iron calorimeter. 

R = 0.65 ±0.05 ±0.08 (Kamiokande) (7.202) 
R = 0.54 ±0.05 ±0.11 (1MB) (7.203) 
R = 0.61 ±0.15 ±0.05 (Soudan2) (7.204) 

The fact that the double ratio R is significantly less than one is an indication in favor of 
neutrino oscillations. 
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The important evidence in favour of neutrino oscillations was obtained by the Super- 
Kamiokande collaboration. These data were first reported at NEUTRIN098 conference 
in Japan, in June 1998. A significant up-down asymmetry of multi-GeV muon events 
was discovered in the Super-Kamiokande experiment. 

For atmospheric neutrinos the distance between production region and detector 
changes from about 20 km for down-going neutrinos (6 = 0, 9 being the zenith angle) up 
to about 13,000 km for up-going neutrinos (9 — ir). In the Super-Kamiokande experi- 
ment for the multi-GeV events the zenith angle 9 can be determined. In fact, charged 
leptons follow the direction of neutrinos (the averaged angle between the charged lepton 
and the neutrino is 15° — 20°). The possible source of the zenith angle dependence of 
neutrino events is the magnetic field of the Earth. However, for neutrinos with energies 
larger than 2-3 GeV, within a few % no ^-dependence of neutrino events is expected. 

The Super-Kamiokande collaboration found a significant zenith angle dependence of 
the multi-GeV muon neutrinos. For the integral up-down asymmetry of multi-GeV 
muon neutrinos (FC and PC) the following value was obtained 

= 0.311 ±0.043 ±0.010 (7.205) 

Here 

A -V^ < 7m > 

where U is the number of up-going neutrinos (cos 9 < —0.2) and D is the number of 
down-going neutrinos ((cos6> > 0.2). No asymmetry of the electron neutrinos was found: 

A e = 0.036 ±0.067 ±0.02 (7.207) 

The Super-Kamiokande data can be described if we assume that there are — > v T 
oscillations. The following best-fit values of the oscillation parameters were found from 
the analysis of FC events 

Am 2 = 3.05 • lO'hV 2 , sin 2 29 = 0.995 (7.208) 



(Xmin = 55-4 at 67 d.o.f.). Let us notice that if we assume that there are no oscillations, 
then in this case x 2 — 177 at 69 d.o.f. From the combined analysis of all data it was 
found 

Am 2 ~ (2 - 6) • 10" 3 e\/ 2 , sin 2 29 > 0.84 (7.209) 

If — > v s oscillations are assumed, at large energies matter effects must be impor- 
tant. From the investigation of the high energy events (PC and upward-going muon 
events, muons being produced by neutrinos in the rock under the detector) the Super- 
Kamiokande collaboration came to the conclusion that — > v s oscillations are dis- 
favoured at 95% C.L. 

The range of oscillation parameters which was obtained from the analysis of the 
atmospheric neutrino data will be investigated in details in long-baseline experiments. 
The results of the first LBL reactor experiment, CHOOZ, were recently published (in 
this experiment the distance between reactors and detector is ~ 1 km). No indication 
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in favour ofthe transitions of v e into other states was found in this experiment. For the 
ratio R of the number of measured and expected events it was found 

R = 1.01 ± 2.8% (stat) ± 2.7% (syst) (7.210) 

These data allow to exclude Am 2 > 7 ■ 10~ 4 eV 2 at sin 2 29 = 1 (90% C.L.). 

In LBL Kam-Land experiment F^'s from reactors at the distance of 150 — 200 km from 
the detector will be detected. Neutrino oscillations V e <-> V x with Am 2 > 10~ 5 eV 2 and 
large values of sin 2 29 will be explored. The BOREXINO collaboration plans to detect 
V e from reactors at the distance of about 800 km from the detector. 

The first LBL accelerator experiment K2K is running now. In this experiment i^'s 
with average energy of 1.4 GeV, produced at KEK accelerator, will be detected in the 
Super-Kamiokande detector (at a the distance of about 250 km). The disappearance 
channel — > and the appearance channel — > u e will be investigated in detail. This 
experiment will be sensitive to Am 2 > 2 • 10~ 3 eV 2 at large sin 2 29. 

The LBL MINOS experiment between Fermilab and Soudan (the distance is of about 
730 km) is under the construction. In this experiment all the possible channels of 
transitions will be investigated in the atmospheric neutrino range of Am 2 . 

The LBL CERN-Gran Sasso experiments (the distance is of about 730 km) ICARUS, 
NOE and others, are under constraction at CERN and Gran Sasso. The direct detection 
of r's from — > v T transition will be one of the major goal of these experiments. 

7.3 LSND experiment 

Some indications in favour of <-> v e oscillations were found in short-baseline LSND 
accelerator experiment. This experiment was done at the Los Alamos linear accelerator 
(with protons of 800 MeV energy). This is a beam-stop experiment: most of 7r + 's in 
the beam, produced by protons, come to a rest in the target and decay (mainly by 
7r + — > ji + Vp); /x + 's also come to a rest in the target and decay by fj, + — > e + v e T7^. Thus, 
the beam-stop target is the source of u^, u e and z/ M (no V e are produced in the decays). 

The large scintillator neutrino detector LSND was located at a distance of about 30 m 
from the neutrino source. In the detector z/ e 's were searched for through the observation 
of the process 

V e +p -> e + + n (7.211) 

Both e + and delayed 2.2 MeV 7's from the capture np — > were detected. 

In the LSND experiment 33.9 ±8.0 events were observed in the interval of e + energies 
30 < E < 60 MeV. Assuming that these events are due to F M — > V e transitions, for the 
transition probability it was found 

Ptpv -> Ve) = (0.31 ± 0.09 ± 0.06) • 10" 3 (7.212) 

From the analysis of LSND data the allowed region in sin 29 - Am 2 plot was obtained. 
If the results of SBL reactor experiments and SBL accelerator experiments on the search 
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for Vp, — > v e transitions are taken into account for the allowed values of the oscillation 
parameters it was found 

0.2 < Am 2 < 2eV 2 2 • 10~ 3 < sin 2 26 < 4 • 10~ 2 (7.213) 

The indications in favour of — > v e oscillations obtained in the LSND experiment will 
be checked by BOONE experiment at Fermilab, scheduled for 2001-2002. 

8 Conclusions 

The problem of neutrino masses and mixing is the central problem of today's neutrino 
physics. More than 40 different experiments all over the world are dedicated to the inves- 
tigation of this problem and many new experiments are in preparation. The investigation 
of the properties of neutrinos is one of the most important direction in the search for a 
new scale in physics. These investigations will be very important for the understanding 
of the origin of tiny neutrino masses and neutrino mixing which, according to the existing 
data, is very different from CKM quark mixing. 

If all existing data will be confirmed by the future experiments it would mean that 
at least four massive neutrinos exist in nature (in order to to provide three independent 
neutrino mass squared differences: Am 2 olar ~ 10 _5 eV 2 (or 10~ 10 eV 2 ), Am 2 tm ~ 10" 3 eV 2 
and Am^gNj-, ~ leV 2 ). From the phenomenological analysis of all existing data it follows 
that in the spectrum of masses of four massive neutrinos there are two close masses 
separated by the "large" one, by the about 1 eV LSND gap. Taking into account big- 
bang nucleosynthesis constraint on the number of neutrinos it can be shown that the 
dominant transition of the solar neutrinos is v e — > z/ stC riie one and the dominant transition 
of the atmospheric neutrinos is v T . 

If the LSND indication in favour of — > v e oscillations will be not confirmed by 
the future experiments, the mixing of three massive neutrinos with mass hierarchy is 
plausible scenario. 

The nature of massive neutrinos (Dirac or Majorana?) can be determined from the 
experiments on the search for neutrinoless double /3- decay. It can be shown that from the 
existing neutrino oscillation data it follows that effective Majorana mass < m > in the 
case of three massive Majorana neutrinos with mass hierarchy is not larger than 10~ 2 eV 
(the present bound is | < m > | ~ 0.5 eV and the sensitivity of the next generation of 
experiments will be | < m > | ~ 0.1 eV). 

The sensitivity | < m > | ~ 10~ 2 eV is very important problem of experiments on the 
search for neutrinoless double /5-decay. 

I would like to express my deep gratitude to R. Bernabei, W.M. Alberico and S. 
Bilenkaia for their great help in preparing these lecture notes. 
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